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ABSTRACT 


The  velocity  of  sound  in  mixtures  of  water  and  steam  has  been  computed 
from  the  slope  of  the  lines  of  equal  entropy  over  a  wide  range  of  ambient 
conditions  up  to  the  critical  point.  The  sound  velocity  is  lower  in  the 
mixture  than  in  either  constituent.  It  changes  dis continuously  on  transition 
to  the  single  phase  in  wet  steam  with  just  a  trace  of  liquid.  The  velocity  is 
about  10  per  cent  less  than  in  dry  steam.  The  velocity  decreases  monotoni- 
cally  with  increasing  water  content.  In  boiling  water  the  velocity  is  very 
low.  The  smallest  value  occurs  at  reduced  pressure  at  the  triple  point 
(with  no  solid  ice  present);  the  velocity  here  is  as  low  as  0,  012  m/s, 

Rankin e- Hug oniot  curves  were  calculated  for  the  mixture.  These  show 
that  a  shock-like  compression  in  very  wet  steam  or  boiling  water  will  con¬ 
dense  all  the  steam,  giving  rise  to  a  large  change  in  volume  with  very  little 
rise  in  pressure.  The  Hugoniot  curve  lies  very  close  to  a  line  of  constant 
entropy  in  this  case.  For  high  quality  steam  the  compression  will  raise 
the  temperature  of  the  vapor  and  evaporate  the  liquid  phase.  This  occurs 
far  from  constant  entropy  conditions,  indicating  a  strong  loss  mechanism. 

In  boiling  water  the  propagation  was  studied -with  pressure  waves  having 
an  amplitude  exceeding  the  minimum  necessary  to  condense  all  the  steam. 
The  velocity  decreases  monotonically  with  decreasing  driving  pressure, 
approaching  the  small  amplitude  acoustic  velocity  when  the  driving 
pressure  approaches  the  minimum  necessary  to  condense  all  the  steam. 

In  a  straight  pipe  of  boiling  water  true  shocks  were  not  propagated. 

The  rise  time  of  the  compression  wave  steadily  increases  as  it  passes 
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through  the  medium.  On  finally  reaching  a  region  of  no  bubbles,  high 
peak  pressures  are  observed.  These  peaks  were  an  order  of  magnitude 
greater  than  the  original  driving  pressure.  This  amplitude  is  apparently 
due  to  the  formation  of  a  water  hammer. 
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PROPAGATION  OF  PRESSURE  WAVES 


IN  A  MIXTURE  OF  WATER  AND  STEAM 

I.  INTRODUCTION 

The  acoustic  properties  of  reactor  media  are  of  interest  in  the  deter¬ 
mination  of  the  propagation  of  shocks  generated  by  a  malfunction.  The  design 
of  enclosures  must  be  based  on  precise  knowledge  of  the  forces  generated, 
and  the  acoustic  impedance  and  attenuation  of  the  medium. 

The  problem  of  a  liquid-gas  mixture  was  treated  in  a  previous  report* 
giving  calculations  and  measurements  of  the  velocity  of  sound  in  a  water-air 

mixture.  The  propagation  of  shocks  through  such  a  mixture  has  been  studied 
2 

in  England  .  This  present  report  is  concerned  with  single  component  two 
phase  liquid-vapor  mixtures,  specifically  water-steam  mixtures,  in  which 
the  total  mass  of  either  phase  will  change  with  the  changing  pressure  of  a 
passing  acoustic  wave.  The  behavior  therefore  differs  from  liquid-gas 
systems,  velocity  of  propagation  being  exceedingly  low  in  the  boiling  fluid. 

II.  COMPUTATION  OF  PRESSURE  PROPAGATION 
IN  A  WATER-STEAM  MIXTURE 

In  a  two  phase  medium  such  as  steam  and  water,  the  propagation 
constants  of  a  pressure  wave  depend  on  the  thermodynamic  variables  of  the 
two  states:  latent  heat,  specific  heat  and  density.  A  wa\e  passing  through  a 
water-air  mixture,  for  example,  has  a  velocity  determined  by  the  average 
compressibility  and  density  of  the  constituents.  As  a  result,  the  medium 

behaves  like  a  fluid  having  a  high  density  like  water  and  a  high  compressibility  ! 

T~ - - - — 

Superscript  numerals  refer  to  bibliography  section. 
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like  air.  The  pressure  wave  compresses  the  air  because  of  its  relatively 
lower  stiffness ,  but  not  the  water.  Also,  the  compression  is  practically 
isothermal  because  of  the  enormous  heat  capacity  of  the  water.  If  the  gas 
is  a  vapor  intermixed  with  its  liquid  phase  at  equilibrium,  pressure  is 
uniquely  related  to  temperature  so  that  if  equilibrium  is  maintained,  a 
change  in  pressure  produces  a  change  in  temperature  and  relative  mass 
and  volume  of  each  phase  through  condensation  and  evapor  ation.  Sound 
waves  of  small  pressure  amplitude  propagate  through  the  vapor -liquid 
medium  isentropically.  Intermediate  amplitude  and  shock  waves  propagate 
nonisentropically  with  dissipation  of  appreciable  amounts  of  energy.  The 
following  discussion  begins  with  small  amplitude  wave  propagation  in  an 
intermixed  vapor -liquid  medium,  in  which  account  is  taken  of  the  change 
of  state  which  occurs  during  the  pressure  cycle. 

A.  Small  Amplitude  Sound  Waves 
1.  Propagation  Velocity 

The  velocity  of  small  amplitude  sound  waves  in  liquid-vapor  mixtures 
is  derived  on  the  assumption  of  instantaneous  equilibrium  conditions  at  all 
points  in  the  system.  If  the  mixing  of  the  two  componert  s  is  intimate  enough, 
equilibrium  will  rapidly  be  attained.  Put  another  way,  for  any  given  particle 
size  of  each  phase  there  will  exist  a  time  in  which  equilibrium  will  be 
established.  Frequency  much  smaller  than  the  reciprocal  of  this  time 
will  be  propagated  under  equilibrium  conditions.  We  are  concerned  here 
with  this  low-frequency  approximation.  The  value  of  this  limit  is  discussed 
in  Appendix  H. 


Report  No. 

.  ARF  4132-12 


In  a  liquid-vapor  equilibrium  system,  pressure  and  temperature 
are  uniquely  related.  For  each  value  of  pressure  there  exists  one  and 
only  one  value  of  temperature,  so  that  a  rise  in  pressure  is  accompanied 
by  a  rise  in  temperature.  This  change  in  conditions  causes  some  evapor¬ 
ation  or  condensation  to  take  place  producing  a  relatively  large  change  in 
volume. 

Now  if  we  take  into  account  the  volume  change  due  to  the  change  in 
pressure,  of  the  vapor,  and  the  change  in  mass  of  vapor  due  to  .'evaporation 
or  condensation,  a  total  rate  of  change  of  volume  (or  its  reciprocal: 
density)  with  pressure  may  be  derived.  In  Appendix  A  we  have  derived 
this  rate  of  change  of  pressure,  p,  with  density,  p,  at  constant  entropy, 

S  ,  analytically.  This  rate  of  change  is  equal  to  the  square  of  the  velocity 

of  sound,  c:  ' 

c2  =  (3p/8p)s  .  (1) 

In  terms  of  specific  volume  of  fluid,  Vf  ;  specific  entropy  of  fluid,  Sf  ; 
the  changes  of  these  quantities  upon  evaporation,  V  and  Sg  ;  and  the 
quality,  q  (relative  mass  of  steam),  the  velocity  is  given  by 


c  =  (q  Ve  +  Vf) 


dS 


L  Se 


d  p 


dV  V  dS,  dV, 
_ e  _e  _ f _ f 

dp  Sg  dp  “  dp 


1-1/2 


(2) 


Let  us  look  at  pressure-density' curves ,  Fig.  1,  and  see  how  the 
slopes  of  these  lines  change  with  given  conditions.  We  observe  that  the 
lines,  though  of  necessity  continuous  through  the  liquid-mixture  and 
mixture-vapor  transitions,  do  not  have  continuous  slopes;  the  lines 
actually  kink  sharply,  having  a  lower  slope  in  the  mixture  than  in  either 
phase.  This  slope  (if  plotted  on  linear  instead  of  logarithmic  paper)  is 
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equal  to  the  square  of  the  sound  velocity.  It  is  seen  that  the  slopes  are 
greatest  in  pure  water  on  the  right  and  least  in  the  mixture  immediately 
adjacent  to  it,  that  is,  in  boiling  water  containing  a  small  quantity  of 
steam.  The  very  low  velocity  of  sound  in  boiling  water  with  low  steam 
content  may  be  understood  by  realizing  that  a  very  small  pressure  change 
condenses  a  large  fraction  of  the  available  vapor  to  give  a  large  percentage 
change  in  the  over-all  density.  We  are,  of  course,  assuming  an  acoustic 
pressure  that  is  so  small  that  the  vapor  content  will  decrease  without 
vanishing  during  the  positive  half  cycles.  Similar  limiting  considerations 
apply  to  the  discontinuous  velocity  change  when  the  vapor -fluid  mixture 
approaches  the  pure  vapor  region.  It  is  of  interest  to  observe  in  Fig.  1 
that  a  pressure  increase  in  high  quality  steam  (relative  mass  of  steam 
large)  brings  the  mixture  closer  to  the  pure  vapor  region.  This  is  due 
to  the  adiabatic  temperature  rise  which  causes  some  of  the  liquid  to 
evaporate.  When  the  mass  of  liquid  is  large,  this  heat  is  absorbed  and 
a  pressure  increase  condenses  some  steam.  The  dividing  line  between 
evaporation  and  condensation  for  a  w&ter-steam  mixture  is  close  to  the 
50  per  cent  concentration  at  all  ambient  pressures.  If  water  droplets 
in  steam  or  vapor  bubbles  in  liquid  have  exceedingly  small  radii  (an 
assumption  necessary  to  extend  the  range  of  the  low-frequency  approxi¬ 
mation),  then  the  effect  of  surface  tension  will  become  noticeable.  In 
that  case  equilibrium  conditions  demand  a  difference  in  pressure  across 
the  curved  interface.  This  differential  has  been  neglected  in  the  present 
treatment,  limiting  the  low-frequency  approximation. 

The  velocity  of  sound  as  given  by  Eq.  (2)  was  computed  for  many 
values  of  pressure  (or  temperature)  and  quality.  The  values  of  specific 
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volume  and  entropy  were  taken  from  Keenan  and  Keyes:  Thermodynamic 

3 

Properties  of  Steam  .  The  derivatives  were  obtained  by  using  a  chord 
approximation,  taking  values  from  alternate  lines  in  the  steam  tables. 

The  errors  introduced  by  this  technique  are  quite  small  as  shown  in 
Appendix  B.  The  variation  of  velocity  with  temperature  and  quality  is 
shown  in  Fig.  2.  A  somewhat  similar  graph  is  shown  as  Fig.  3  with  a 
change  of  variable  for  the  abscissa  which  replaces  temperature  with  the 
logarithm  of  the  pressure.  With  this  change  of  variable  the  sound 
velocity  approximates  a  linear  function  of  the  independent  variable,  the 
logarithm  of  pressure,  at  very  low  steam  concentrations.  Individual 
curves  correspond  to  specified  values  of  the  quality  of  steam,  q.  In 
Fig.  4  velocity  has  been  plotted  as  a  function  of  quality  or  phase  concen¬ 
tration  at  selected  values  of  pressure,  p  (or  the  corresponding  values 
of  temperature,  T).  A  monotonic  increase  in  velocity  with  concentration 

for  all  pressure  is  exhibited.  The  rate  of  increase  of  velocity  with 

-  2 

quality,  q,  is  largest  for  a  pressure  in  the  region  of  100  lb  in.  while 
it  is  smaller  for  pressures  both  below  and  above  this  value.  As  the 
critical  point  is  approached,  the  energy  required  for  a  phase  change 
vanishes,  so  that  concentration  becomes  meaningless  and  the  sound 
velocity -quality  curve  approaches  a  horizontal  line.  In  Fig.  5  the  same 
curves  are  plotted  on  a  logarithmic  scale.  Here  the  variations  at  the 
lower  concentrations  are  emphasized.  We  note  that  in  water  containing 
just  a  few  steam  bubbles,  a  very  low  sound  velocity  compared  to  the 
velocity  in  pure  water  is  derived.  The  change  to  the  high  velocity  in 
pure  water  is  discontinuous. 
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The  computations  have  been  extended  to  the  condition  of  pure  dry- 
steam,  and  the  results  are  presented  in  graphical  form  in  Fig.  6.  It  is 
observed  that  at  the  same  temperature  and  pressure  pure  dry  steam  has 
a  higher  velocity  than  wet  steam,  even  when  the  water  content  tends  to 
zero;  that  is  to  say,  we  compute  a  discontinuous  jump  in  the  velocity  at 
this  point  also.  This  is  clearly  shown  in  Fig.  6,  in  which  velocity  is 
plotted  as  a  function  of  pressure  for  constant  entropy.  Lines  of  constant 
quality  and  lines  of  constant  temperature  are  also  shown.  At  low 
pressures  and  high  temperatures,  velocity  is  reasonably  independent  of 
pressure  as  one  would  expect  for  a  perfect  gas.  At  higher  pressures 
3-rid  lower  temperatures  these  lines  tend  to  curve  down.  It  is  interesting 
to  note  here  that  if  instead  of  plotting  velocity  as  a  function  of  temper¬ 
ature,  we  plot  it  as  a  function  of  constant  enthalpy,  as  in  Fig.  7,  the 
velocity  is  almost  independent  of  pressure,  even  up  to  very  high  pressures, 
provided  condensation  does  not  take  place;  see  also  Fig.  8.  When  con¬ 
densation  does  take  place  we  can  observe  discontinuous  velocity  changes 
when  the  slightest  trace  of  condensation  occurs.  The  region  of  very  high 
pressures  and  low  enthalpy,  in  which  the  lines  are  seen  not  to  be  per¬ 
fectly  horizontal,  is  rather  uncertain  because  the  points  tabulated  in  the 
steam  tables  are  too  far  apart  for  a  determination  of  the  differential 
parameters  as  required  by  Eq.  (1).  The  difference  between  the  water- 
steam  mixture  and  water-air  mixture,  covered  in  our  previous  study^, 
is  brought  out  in  Fig.  6,  where  we  see  that  for  constant  entropy  the  velocity 
of  sound  rises  very  slightly  with  pressure  at  high  entropy,  more  slowly 
in  fact  than  in  a  single  phase  medium,  and  at  low  entropy  or  low  quality, 
the  velocity  actually  decreases  with  increasing  pressure. 
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2. 


Comparison  with  Data  from  the  Literature 


There  have  been  relatively  few  measurements  of  the  velocity  of  sound 

4 

in  steam.  Velocity  was  measured  by  a  resonance  technique  by  Masson  in 

C 

1857;  Jaeger  in  1889  used  a  Kundt's  tube  technique;  and  Treitz  reported  on 

the  velocity  of  sound  in  steam  in  1903.  Comparison  of  the  observations  with 

our  computation  is  summarized  in  Table  L  It  is  seen  that  the  results  of 

7 

Jaeger  and  of  Treitz  are  at  variance  with  those  of  Shilling  (1924)  and  those 

Q 

of  Billhartz  and  Bishop  (1936).  Shilling  was  very  careful  to  remove  liquid 

water  droplets  and  also  obtained  his  results  at  the  end  point  or  extrapolation 

of  a  series  of  measurements  in  superheated  steam. 

It  is  very  tempting  to  assume  that  the  steam  used  by  Jaeger  and  by 

Treitz  contained  water  droplets.  For  the  measurements  of  saturated  steam 

(at  100° C  and  below)  an  assumption  of  about  14  per  cent  by  weight  of  water 

would  bring  these  observations  in  line  with  computations.  This.,  of  course, 

does  not  explain  Treitz'  results  at  110,  120,  130°C  at  atmospheric  pressure. 

To  explain  these  we  postulate  that  even  in  slightly  superheated  steam 

liquid  droplets  will  exist  .  Surface  tension  provides  an  excess  pressure 

inside  the  droplet  so  that  it  is  in  equilibrium  with  the  vapor  at  the  same 

temperature.  This  leads  to  a  stable  droplet  of  radius  r  =  2  T/(p^  -  pa), 

where  T  is  the  surface  tension  and  (p^  -  pa)  is  the  difference  between 

ambient  pressure  and  boiling  pressure  at  that  temperature.  In  the  case 

-4 

of  Treitz'  work  droplet  radii  of  the  order  of  10  cm  are  calculated.  These 
may  well  escape  observation,  being  no  larger  than  a  wavelength  of  light. 

Changes  in  diameter  will  of  course  have  an  effect  on  the  equilibrium 
conditions  because  of  the  change  in  surface  energy.  The  smaller  the 
radius  the  larger  the  error  introduced  by  neglecting  this  effect.  For 
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droplets  with  radii  greater  than  2  T/Lp  (T  =  surface  tension,  dyne  cm-2, 

-  1  * 

L,  =  latent  heat,  dyne  cm  gm  ,  p  -  density,  gm  cm"  )  the  surface  tension 
effect  is  smaller  than  the  latent  heat  effect.  Even  for  the  droplets  of  10~^ 
cm  radii  postulated  above,  surface  tension  has  a  negligible  effect  on  the 
sound  velocity  and  is  justifiably  neglected. 

Evidence  for  the  existence  of  liquid  droplets  in  superheated  vapor 
is  found  in  the  effectiveness  of  and  necessity  for  the  special  filters 
Shilling  used  to  remove  droplets  during  his  measurements,  all  of  which 
were  at  atmospheric  pressure  and  above  100* C.  It  is  also  worth  noting 
that  the  above  relation  for  the  radius,  r  =  2  T/(p,  -  p  ),  yields  zero  radius 
above  the  critical  temperature  when  the  surface  tension  vanishes.  This 
means  that  there  are  no  droplets  in  "permanent"  gases  (i.  e.  ,  gases  above 
the  critical  temperature  which  cannot  be  liquefied  by  simple  compression). 

The  formula  also  yields  an  infinite  radius  for  droplets  in  saturated 
steam.  Frenkel10  postulates  that  droplets  in  saturated  steam  are  hotter 
than  the  average  temperature  of  the  vapor  and  are  surrounded  by  a  thin 
layer  of  superheated  vapor.  As  they  lose  heat  by  conduction,  they  grow 
in  size  and  will  eventually  fall. 

Droplets  in  superheated  steam  will  remain  the  equilibrium  size, 
larger  ones  evaporating  and  smaller  ones  gathering  more  condensation. 
Ultimately  collision  between  droplets  and  between  droplets  and  walls 
causes  the  larger  agglomerate  to  evaporate  down  to  the  stable  size,  thus 
slowly  eliminating  the  liquid  phase. 
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Table  1 


THE  VELOCITY  OF  SOUND  IN  STEAM 


Observer 

Temp. 

°C 

c  reported 
m/ s 

c  calculated 
wet*  dry** 

q  -*•-  0 

Mas  son ^ 

0 

402 

Jaeger^ 

93 

402 

434 

469 

Jaeger 

96 

410 

437 

472 

5 

Treitz 

100 

404.  3 

439 

475 

Treitz 

110 

413.  0 

444 

478 

Treitz 

120 

417.  5 

449 

482 

Treitz 

130 

424.  4 

455 

488 

Shilling  ^ 

100 

471.  5 

439 

473.  1 

Shilling 

300 

593.  2 

❖❖❖ 

586 

Shilling 

600 

727.  8 

❖❖❖ 

714 

Shilling 

800 

795.  3 

5!'  * '!' 

787 

Billhartz  ^ 

100 

471.  3 

439 

473.  1 

and  Bishop 

According  to  Eq.  (2). 

❖❖  2 

According  to  c  =  \  p  V  ;  values  of  y  and  V  were  taken 

cfaam  f  oKl  o  e  i  S  S 


*  >:=  ❖ 


from  steam  tables 
Liquid  not  stable  at  atmospheric  pressure  at  this  temperature. 
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B.  An  Empirical  Expression  for  the  Pressure-Volume 
Relation  at  Constant  Entropy 

An  expression  has  been  obtained  which  gives  the  specific  volume  V 
of  a  water-steam  mixture  in  terms  of  specific  entropy,  S,  and  pressure,  p: 

V  =  Aq  (S  +  0.  24)  p* 1  +  Bq  (S  -  0.  76)  p"°*  7.  (3) 

Values  of  Aq  and  Bq  are  given  for  several  systems  of  units  in  Table  2. 
The  expression  is  an  approximate  fit  of  the  measured  values  of  entropy  and 
pressure.  Accuracy  of  the  parameters  Aq  and  Bq  is  better  than  10  per 

_  O 

cent  in  the  range  0.  1  <  p  <  3200  lb  in,  (ice  point  to  critical  point),  and 

is  better  than  4  per  cent  over  the  most  important  part  of  the  range  between 

-2 

15  and  2500  lb  in.  .  However,  when  S  is  less  than  about  0.  5,  the  value 
for  V  is  the  small  difference  between  large  quantities,  and  the  expression 
is  not  reliable.  Further  work  is  needed  to  obtain  a  good  expression  for 
small  entropies. 

A  functional  relation  of  this  kind  greatly  facilitates  differentiation  to 
obtain  an  approximate  value  for  sound  velocity,  and  integration  under  an 
adiabatic  line  to  obtain  work  done  if  no  heat  is  lost.  The  form  of  the 
function  is  relatively  easy  to  handle. 

The  velocity  of  sound  is 

c  =  v  (-9p/av)s1/2  , 

r  -2  1 71  -1/2 

=  V  Aq  (S  +  0.24)p  S0.7  Bq  (S  -0.  76)p" i*  1  .(4) 
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The  work  done  going  from  condition  Vj  to  V2  at  constant  entropy 
is 


Pi 

=  A  ^(S  +  O.  24)  log  (pj/p^J  +  B  ^S-O.  76)(p1°*  3-p2°‘  3)>  (5) 

where  A,  =  A  and  B,  =  2.  33  B  in  a  consistent  system  of  units.  For 
1  o  1  o 

other  systems  see  Table  2« 

The  method  of  obtaining  Equation  (3),  the  expression  for  volume  in 
terms  of  pressure,  is  detailed  in  Appendix  D* 

C.  Rankine-Hugoniot  Curves  for  Pressures  up  to 
Complete  Condensation  or  Evaporation 

The  change  of  specific  volume  V  and  specific  internal  energy  E 
across  a  shock  front  are  given  in  terms  of  the  pressure  p  by  the  well- 
known  Rankine-Hugoniot  relationship, 

E2  -  Ej  =  (1/2)  (px  +  p2)  (Vl  -  V2),  (6) 

where  the  subscripts  1  and  2  refer  to  conditions  before  and  after  the  shock. 
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Table  2 


PARAMETERS  FOR  THE  PRESSURE^  VOLUME  RELATION 


V  =  Ao  (S  +  0.  24)  p"1  +  BQ 

(S  -  0, 

76)  p"0<  7 

Units  of  Parameters 

Specific  Energy 

P  V 

A 

o 

Bo 

A! 

B1 

lb  in.  3  ft ^  lb"* 

lb  in.  ~  ^ft3lb"  * 

113 

67.  5 

A 

A 

157 

lb  in.  "  3  ft3  lb"  * 

Btu  lb'1 

113 

67.  5 

U 

610 

847 

,  -2  2-1 
dyne  cm  cm  gm 

-1 

erg  gm  4. 

9xl08 

1.  OxlO7 

A 

o 

2.  4xl07 

newton  m  3  m3/k  g 

joule  ,  kgm"1  4. 

9xl04 

2.  OxlO3 

Ao 

4.  8x10  3 

atmospheres  cm^gm"* 

cal  gm  1 

480 

640 

2.  OxlO4 

6.  OxlO4 

Note:  For  the.;specific- entropy  S,  1  Btu/.*F  lb  =  1  cal/°C  gm; 
hence  Sq  =  0.  24Btu/°Flb 

=  0.  24  cal/°C  gm. 
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In  the  case  of  a  water-steam  mixture  we  get  the  additional  four 
relationships 

V1  =  ql  Vel  +  Vfl  ’  V2  =  q2  Ve2  +  Vf2  ’ 

and 

E1  =  ql  Eel  +  Efl  ’  E2  =  q2  Ee2  +  Ef2  * 
for  conditions  1  and  2. 

For  an  equilibrium  process  (that  is,  after  completion  of  relaxations 
near  the  wave  front)  E^  and  q^are  eliminated  between  these  equations. 

We  obtain  then 

(Pl  +  P2>  V1  +  2G  <E1  -  E12  +  Ee2  Vf2/Ve2> 

V2  =  (Pl  +  p2>  +  2G  (Ee2/Ve2)  1 

where  the  constant  G  is  unity  for  a  consistent  system  of  units.  In  our  case 

_  2  3-1  .  - 1 

where  steam  tables  give  p  in  lb  in  ,  V  in  ft  lb  ,  and  E  in  Btu  lb  the 
factor  G  is  5.  4046. 

Two  Rankine -Hugoniot  curves  are  shown  in  Fig.  9,  where  has 

been  plotted  against  p2-  For  the  first  case,  curve  (a),  p1  has  been  chosen 
2  3-1 

as  14.  7  lb  in.  ”  and  Vj  as  1  ft  lb  .  For  the  second,  curve  (b),  p^  is 

14.  7  lb  in.  "2  and  Vj  is  10  ft3  lb-1.  We  note  that  the  quantity  Ej  is 

determined  once  p^  and  have  been  specified,  and  that  the  quantities 

E  ,,  V,-,,  and  V  0  are  constants  of  the  fluid-gas  mixture.  The  curve 
e2  f2  e2 

constructed  of  short  dashes  gives  the  pairs  of  values  for  pressure  and 
specific  volume  at  constant  entropy. 

-  2 

We  observe  that  for  curve  (a),  when  a  pressure  of  29  lb  in.  or 
about  2  atmospheres  is  reached,  all  steam  has  been  condensed  to  water 
and  the  volume  has  been  reduced  to  l/60th  of  its  original  value.  This 
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would  require  a  large  amount  of  energy.  We  observe  further  that  if 
we  draw  an  adiabat  through  the  final  point,  it  almost  exactly  coincides 
with  the  Rankine "Hugoniot  curve.  A  similar  curve  could  be  plotted 
starting  with  a  specific  volume  of  0.  1  ft  /lb  in  which  case  the  pressure 
reached  for  complete  condensation  does  not  exceed  16  lb  in.  That 
is,  a  sixfold  volume  reduction  would  be  observed  with  less  than  a  10  per 
cent  pressure  increase.  In  general,  in  bubbly  water  (less  than  50  per 
cent  steam  by  volume)  a  very  small  pressure  increase  will  condense  all 
the  steam.  The  energy  required  to  do  this  (since  pressure  is  essentially 
constant)  is  of  course  simply  p  V,  where  V  is  the  total  volume  of  steam. 
Thus  any  explosion  with  less  than  this  amount  of  energy  could  never 
produce  a  significant  pressure  increase  in  an  equilibrium  process.  The 
value  of  the  minimum  pressure  needed  to  condense  all  the  steam  in  a 
mixture  is  discussed  in  Section  II- D- 3  and  in  Appendix  E. 

If  we  were  to  consider,  for  example,  an  initial  condition  of  250  ft^ 

of  water-steam  mixture  containing  20  per  cent  steam  by  volume  at  a 

-  2 

pressure  of  700  lb  in.  ,  the  pressure  increase  required  to  condense  all 
the  steam  adiabatically  may  be  arrived  at  as  follows: 

3 

volume  of  water  200  ft 

3 

volume  of  steam  50  ft 

entropy  1.  425  Btu/lb°  F 

Interpolation  in  the  steam  tables  show  that  this  entropy  corresponds  to 
water  without  steam  at  725  lb  in.  pressure. 

Thus  in  this  case  there  is  a  pressure  rise  of  about  3.  5  per  cent.  To 
obtain  the  energy  required  to  perform  this  compression  it  is  a  fair  approxi¬ 
mation  to  consider  the  work  done  in  compressing  this  steam  as  "p  V,  where 
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p  is  the  average  pressure,  or  712  lb  in. 


-  2  3 

W  =  712  x  50  lb  in,  ft  , 

=  6590  Btu  , 

=  1,  65  x  106  cal, 

i  j  mixture  in  this  example  is 

This  shows  that  if  the  energy  released  in  the  mixtur 

less  than  ahont  ,  6  million  eateries,  there  witt  he  very  little  eventual 

pressure  build-up  at  the  wall  of  the  containing  vessel.  On  the  o 

4.1  •  of  this  pressure  amplification 

if  the  energy  release  is  great  y  m 

may  occur,  as  shown  in  Section  II- D, 

The  other  case,  illustrated  in  curve  <b>  of  Fig.  9,  starts  with  we, 

steam  of  specific  volume  10  ft3  lb'  1  at  one  atmosphere  of  pressure.  In 
this  case  the  shock  is  seen  to  evaporate  water  until,  at  a  pressure  of 
about  680  lb  in.  pure  steam  remains.  In  this  case  the  Rankine-Hugomot 

curves  and  the  adiabat  are  separated  leaving  the  steam  in  a  state  of  higher 
quality,  at  a  higher  pressure  in  a  given  volume,  and  with  higher  sp 
energy.  Energy  is  obviously  required  to  increase  the  energy  m  this 
manner  and  the  only  possibie  source  is  the  shock  wave  itself.  The  energy 
of  the  shock  must  therefore  be  attenuated  in  traveling  through  the  medium. 

„„  th.  energv  under  the  adiabat,  Eq.  (5)  of  Section 
In  order  to  compute  the  energy 

.  -2  14  7  1b  in  ^  s  =  1,  433 , 

H.B  is  useful.  For  p2  =  680  lb  in.  ,  Px  14.  7  lb  in. 

680 


-141  1  t>eu  4. 

w  =  610  (1.  433  +  o,  24)  loge 
847  (i.  433  -  0,  76)  (680 


.0,  3 


14.  7°'  3), 


=  6673  Btu /lb. 


-15- 


Report  No. 
ARF  4132-12 


D. 


Propagation  of  Large  Amplitude  Pressure  Waves 


1.  General  Description  of  Events 

The  preceding  section  described  the  behavior  of  shocks  of  small 
amplitude,  so  that  the  mixed  medium  remained  a  mixture  after  the 
event.  In  boiling  liquid  in  which  the  original  density  is  within  25  per  cent 
of  that  of  the  liquid  the  pressure  required  to  condense  all  the  steam  was 
shown  to  be  very  small  (Fig.  9;  see  also  Appendix  E).  We  now  consider 
larger  pressures  which  condense  all  the  steam  to  a  single  liquid  state. 
These  were  also  studied  experimentally  (see  Section  III). 

The  application  of  a  pressure  pulse  to  the  boiling  liquid  causes  the 
steam  bubbles  near  the  surface  to  collapse,  and  a  wave  of  bubble  collapse 
travels  away  from  the  point  of  application  of  the  pressure  until  all  the 
vapor  is  condensed.  The  events  are  shown  schematically  in  Figs.  10,  11 
and  12. for  an  arrangement  in  which  a  shock  wave  is  applied  from  above 
to  the  boiling  water.  Figure  10A  shows  the  appearance  of  the  liquid- 
vapor  mixture  extending  from  the  surface  down  to  some  point  m  before 
the  pulse  is  applied.  The  region  below  m  is  filled  with  pure  liquid. 
Some  time  after  a  pulse  has  been  applied  (Fig,  10B)  the  mixture  in  the 
region  e  -  j  has  been  condensed  to  pure  liquid,  and  the  surface  has  been 
displaced  from  d  to  e.  By  the  time  the  front  has  reached  m,  the  end 
surface  is  at  point  h,  below  the  original  level  at  d.  First,  we  consider 
the  events  which  would  be  observed  with  a  pressure  transducer  at 
position  c  in  Fig.  10  just  above  the  liquid  surface.  At  a  time  t^  =  Lbc/cl 
after  release  of  a  pressure  wave  at  b,  an  abrupt  pressure  increase  will 
be  observed  (Fig.  11).  Here  c^  is  the  velocity  of  the  shockwave  in 
the  gas  above  the  boiling  water,  and  Li^  is  the  distance  from  b  to  c. 
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At  time  =  t1  -I-  L  ,  (i-  +  — )  ,  a  reflection  from  the  surface  produces 
C*  X  CCl  C  j  C  ^ 

another  abrupt  pressure  increase.  Velocities  c^  and  c ^  for  the  incident 
and  the  reflected  shock  waves  differ  only  slightly  from  the  velocity  of 
sound  re  if  small  amplitude  shock  waves  are  used.  For  greater  ease 
of  reading  we  write  c  for  an  average  transport  velocity,  so  that 
'2=tl  +  2  Lcd/c  '  A  small  reflection  from  expansion  around  the  fixtures 
at  the  pressure  release  point  b  is  next  observed  at  time  t^  =  ^  /c. 

While  these  events  have  been  occurring,  a  rarefaction  due  to  sudden  release 
of  pressure  in  region  a  has  been  traveling  up  the  driving  tube.  This 
rarefaction  undergoes  a  reflection  at  the  upper  end  of  the  tube  and  appears 
at  position  c  at  time  t ^  =  (2  )/c.  The  pressure  pulse  which 

passed  point  c  at  time  t^  traveling  upward  is  also  reflected  at  the  top, 
and  is  observed  again  at  c  at  time  t^  =  t^  +  2  L  ^/c*  This  sequence 
then  repeats  i  tself  with  slightly  reduced  amplitude  on  each  succeeding 
pass.  Several  cycles  of  the  event  are  sketched  in  the  top  part  of  Fig*  12„ 
Pres  sure  “time  functions  in  the  liquid  at  the  points  i,  j,  k,  o  in 
Fig,  10  are  also  shown  in  Fig,  1L  At  point  i  immediately  below  the 
surface  the  pressure  rises  to  the  full  amplitude  of  the  reflected  shock  wave 
observed  above  the  surface.  This  is  easily  explained  by  the  boundary 
conditions  on  pressure  and  particle  velocity  at  the  interface  between  air 
and  the  dense  water  medium.  The  time  of  transit  through  the  liquid  to 
succeeding  positions  down  the  tube  is  evident.  The  attenuation  of  the  wave 
also  shows  up  as  a  reduction  of  the  wave  front  slope.  Finally*  at  time 
^6  ~  V  *  -^dm^ Cf  t^le  ^nter^ace  between  liquid  and  mixture  has  arrived 
at  point  m,  and  this  instant  a  large  pressure  pulse  is  observed  everywhere 
in  the  liquid.  The  velocity  of  this  pulse  in  the  mixture  is  denoted  by  c^; 
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h  is  about  half  way  between  t^  and  t^.  The  maximum  amplitude  of 
this  pulse  is  much  larger  near  the  bottom  of  the  mixture  column  than 
it  is  at  the  surface. 

The  details  of  this  pressure  pulse  are  shown  in  Figs.  13  and  14. 
Assumed  here  are  a  column  of  mixture  of  length  a  supported  by  a 
column  of  pure  liquid  of  length  b,  which  in  turn  is  terminated  by  a 
rigid  wall.  The  two  cases  of  Figs.  13  and  14  show  the  slightly  different 
behavior  observed  for  a  greater  than  and  less  than  b. 

Just  before  the  last  few  bubbles  collapse,  we  essentially  see  a 
column  of  liquid  of  length  a,  having  a  velocity  v,  approaching  a 
stationary  column  of  length  b.  A  very  short  time  t  after  impact,  the 
particles  in  the  impact  region  will  travel  with  an  intermediate  velocity. 
Conservation  of  momentum  (or  symmetry)  considerations  dictate  that 
with  reference  to  this  transition  region  the  incoming  velocities  must 
be  equal,  i.  e.  ,  v/  2.  Thus  a  wave  with  particle  velocity  v/  2  travels 
in  both  direction  from  the  point  of  impact  (Figs.  13A  and  14A).  The 
pressure  increase  in  this  region,  as  shown  in  Figs.  12B  and  13B,  is  equal 
to  p/  2  =  v  c^/ 2,  and  c-  £  being  density  and  sound  velocity  for  the 
pure  water  without  bubbles.  (For  proof  of  this  statement  see  Section 
III-D-3  below.  )  For  the  case  b  less  than  a  as  in  Fig.  13,  the  pulse 
is  first  reflected  at  the  rigid  wall,  doubling  the  pressure  as  shown  at  D. 
This  double  step  pulse  continues  to  travel  to  the  free  surface,  where  each 
step  is  reflected  as  shown  at  E-F,  G-H  to  return  down  as  a  double  step 
pulse  with  negative  particle  velocity  as  at  I  and  pressure  reduced  to  zero. 
The  slight  decay  of  pressure  and  velocity  after  reflection  at  the  rigid  end 
takes  into  account  the  fact  that  the  boundary  is  a  plate  having  large  rigidity 
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but  a  finite  mass.  The  applied  pressure  consequently  starts  it  moving, 
giving  finite  velocities  at  this  end  in  E,  G,  and  I.  The  pulse  traveling 
down  again  would  become  a  negative  pressure  pulse  on  reflection  as  shown 
dotted  at  K-P.  However,  the  liquid  cannot  support  a  negative  pressure, 
especially  near  its  boiling  point,  and  no  negative  pulse  is  observed. 

If  the  length  a  of  the  moving  column  is  less  than  that  of  the  fixed 
column  b  then  events  differ  in  detail  from  those  just  described.  As:  seen 
in  Fig.  14,  the  pulse  traveling  up  toward  the  free  surface  is  reflected 
first  and  a  complete  pulse  with  a  positive  front  and  negative  return  to 
zero  travels  down  toward  the  rigid  boundary.  There  it  is  reflected  and 
gives  rise  to  a  pressure  doubling  at  this  point  for  a  brief  period  until  the 
rear  of  the  pulse  has  also  been  reflected  (F,  H,  J).  This  pulse,  after 
further  reflection  at  the  free  surface,  would  give  rise  to  a  negative  pressure. 
The  consequence  of  this  negative  pulse  will  be  discussed  after  the  pressure- 
time  diagrams  have  been  considered. 

The  pressure-time  diagrams,  Fig.  15,  are  almost  self  explanatory. 
They  were  obtained  by  considering  the  pressure  variation  from  diagram 
to  diagram  in  Figs.  13  and  14,  If  the  moving  column  is  shorter  than  the 
fixed  column,  then  two  distinct  pulses  will  be  observed  (D,E)  except  for 
a  very  small  region  (F,  G)  near  the  rigid  bottom.  The  separation  vanishes 
at  a  distance  a,  equal  to  the  length  of  the  incoming  liquid  column.  Near 
the  surface  the  reflected  wave  tends  to  cancel  the  original  wave  in  part  „ 
because  of  the  relatively  slow  rise  time  of  the  pressure  pulse;  hence  the 
small  amplitude  in  Fig.  15D.  In  the  region  near  the  rigid  bottom  a  double 
step  (Fig.  15F)  is  observed,  the  length  of  the  halfway  plateau  being 
proportional  to  the  distance  from  the  bottom,  vanishing  completely  at 
the  bottom  (se  e  Fig.  15G). 
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When  the  fixed  column  is  longer  than  the  incoming  moving  column, 
it  is  difficult  to  move  out  of  the  region  of  the  double  step  (Fig.  15  B) 
without  moving  into  the  region  of  low  amplitude  (Fig.  15C)  near  the 
surface.  If  the  fixed  column  length  b  is  short  enough,  twD  pulses 
will  not  resolve  at  all. 

The  negative  pulse  previously  considered  is  not  observed  because 
the  liquid  starts  to  boil  again  and  the  pulse  is  transmitted  as  a  large 
reduction  in  density  and  very  small  reduction  in  pressure.  This  "density" 
pulse  travels  through  the  water  just  like  a  pressure  pulse  and  is 
reflected  at  the  surfaces.  The  velocity,  however,  is  governed  by  water  - 
steam  mixture  properties,  so  that  the  pressure  remains  at  ambient  until 
the  positive  reflection  of  the  pulse  in  the  free  surface  has  traveled  all 
the  way  down  to  the  rigid  bottom  and  there,  on  collapsing  the  last  bubble, 
once  again  gives  rise  to  a  large  pressure  as  before.  A  whole  series  of 
large  amplitude  pulses  are  thus  observed  with  long  intervals  of  constant 
ambient  pressure.  The  intervals  get  progressively  shorter  as  the 
amplitude  of  the  pulse  is  attenuated.  This  is  due  to  the  fact  that  the 
negative  pulse  produces  a  smaller  density  reduction  at  each  reflection. 
Experimental  evidence  is  given  in  Section  III-C,  Figs.  2  3,  24. 

2.  Propagation  Velocity  of  the  Compression  Wave 

To  calculate  the  velocity  u  with  which  the  interface  between  the 
original  mixture  and  the  condensed  liquid  above  it  advances  into  the 
mixture,  we  use  the  conservation  of  mass  and  momentum  across  this 
pressure  front.  Observation  has  shown  this  interface  to  be  a  rather  ill- 
defined  region  which  continues  to  broaden  as  it  progresses  through  the 
pipe.  However,  as  far  as  points  outside  the  region  are  concerned  the 
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conservation  laws  hold.  There  is,  of  course,  the  difficulty  of  defining 
the  velocity  of  the  interface  because  of  its  finite,  increasing  thickness. 
However,  the  velocity  u^  of  the  condensed  fluid  behind  it  is  a  perfectly 
well-defined  velocity  at  any  instant,  if  we  consider  a  point  sufficiently 
far  removed  from  the  interface  so  that  this  velocity  is  constant.  Now  the 
velocity  of  advance  of  an  ideal  interface  is  precisely  related  to  the 
velocity  of  the  water  behind;  we  shall  now  compute  these  two  velocities. 

We  define  a  velocity  u2  =  Uj  «.  uf  as  the  relative  velocity  of  the 
interface  with  respect  to  the  condensed  fluid  behind  it;  then  the  two 
conservation  laws  become,  in  terms  of  the  densities  pj,  p2  and  the 
pressures  Pj,  p2  on  the  two  sides  of  the  interface 

Pl  Uj  =  p2  u2  ,  (conservation  of  mass)  (8a) 

Z  2 

p  1  u.  1  "P2u2  =  •  (conservation  of  momentum)  (8b) 

For  convenience  we  now  define  a  steam  volume  concentration,  a  ,  as  the 
volume  of  steam  per  unit  volume  of  mixture.  It  is  related  to  the  mass 
concentration,  q,  used  earlier  by 

(1-1/  Of)  =  (1-1/  q)  Vf/Vg  =  (1-1/  q)  pg/pf  . 

The  density  of  the  mixture  p^  is: 

p!=Q!  Pg+(l-a)  pf  =  (l-a)  pf  if  a<  0.9,  (9) 

since  at  atmospheric  pressure  p^  =  1600  p^.  Thus 

P2  =  pf  =  p  j/(  1-  a)  •  (1°) 
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(Ha) 

(lib) 

(He) 


It  is  of  interest  that  the  expression  for  the  velocity  has  the  same  form 
as  the  velocity  of  sound  in  a  water-air  mixture  ,  except  that  the  ambient 
pressure  p  is  now  replaced  with  the  driving  pressure  (p^  -  Pj). 

Some  caution  must  be  exercised  in  the  use  of  the  formula  for  velocity, 
since  the  derivation  assumed  complete  condensation.  There  exists  a 
pressure  differential  below  which  complete  condensation  does  not  take 
place.  In  Appendix  E  this  minimum  pressure  is  shown  to  be 


<P2  ■  Pl>mm  =  aO-a)'1  vf  Se  Vg‘  1  (dSj /*pf 1  .  (12) 

For  a  =  0.  5  at  100*0  (212’F),  (p2  -  p^^  =  0.  18  psi. 

The  value  of  (p^  -  Pi)m^n  is  plotted  in  Fig.  16. 

It  is  also  shown  in  Appendix  E  that  when  the  pressure  differential  only 

barely  condenses  all  the  water,  i.  e.  ,  when  p  -  p.  =  (p.  -  p_)  .  ,  the 

^  i  c  JL  •  mxn 

velocity  (u.)  is  very  nearly  equal  to  the  small  amplitude  sound  wave  speed 
*  mxn  * 

discussed  in  the  first  part  of  this  report.  Therefore  we  see  that  above  this 
minimum  driving  pressure,  propagation  velocity  decreases  with  the  de¬ 
creasing  driving  pressure  until  this  driving  pressure  is  so  small  that  not 
all  the  steam  is  condensed.  When  this  happens  a  further  reduction  of 
driving  pressure  has  little  effect  on  the  propagation  velocity. 
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3. 


pressure  Buildup  in  a  Water  Hammer 


Consider  two  columns  of  water  approaching  each  other  with 
relative  velocity  u^  in  an  evacuated  pipe.  Since  both  columns  are 
alike  the  velocity  of  the  first  particles  to  touch  after  collision  will  be 
the  mean  of  the  velocities  of  the  two  columns,  i.  e.  ,  u^/ 2.  This  will 

be  their  velocity  with  respect  to  each  column.  Each  column  therefore 
sees  a  pressure  wave  advancing  from  the  point  of  collision  with  a 
particle  velocity  uf/  2.  If  this  wave  travels  along  the  column  with 
velocity  c,  a  length  of  column  AL  =  c  A  t  will  have  stopped  after  the 
time  A  t.  This  length  has  a  mass  pAc  At,  ;The  change  of  momentum 
is  p  A  c  A  t  (u^  -  1/2  u^.),  and  the  time  rate  of  change  of  momentum 
becomes  p  A  c  u^  /  2.  This  is  equal  to  the  total  force  F  =  PA.  Hence 
the  peak  pressure  P  =  p  c  u^/2. 

We  may  further  introduce  the  concept  of  pressure  amplification 
P/ (p^  -  p^)»  when  a  driving  pressure  p^  -  p^  produces  in  the  fluid  a 
velocity  u^  which  at  the  end  of  the  bubble  region  produces  the  water 
hammer  pulse  P. 

Thus 

P  =  Pf  cf  uf / 2  =  (1/ 2)  pf  cfy/(p2  -  p j j Of/  (1  -  Of)  p^»  (13) 

And  the  pressure  amplification  becomes 

r  n1/2 

P/  (p2  -  p£)  =  (1/  2)  c£  |_p£  a!  (1  -  a)  (p2  -  Pj)J  (14) 


=  75]/  Oil  [d  -  a)  (pb2  -  PbiO 


=  287 y  a/ [(1  -  a)  (Pp2-  ppl)]  . 

where  p^  is  the  pressure  measured  in  bars  and  pp  is  the  pressure  measured 
in  psi. 
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Consider  a  collapsing  bubble.  The  radial  velocity  of  the  water 
surrounding  it  endows  this  shell  of  water  with  kinetic  energy.  When 
the  bubble  radius  finally  becomes  zero,  the  radial  velocity  and  the  kinetic 
energy  associated  with  it  vanish.  This  energy  is  converted  to  heat.  To 
calculate  the  energy  dissipated  in  this  manner  we  calculate  the  kinetic 
energy  of  the  mass  of  water  behind  the  "shock"  front  in  our  case  and 
compare  it  with  the  work  done  on  it  by  the  force  accelerating  it. 

Now  from  the  velocity  u^  and  mass  M,  we  may  compute  the  kinetic 
energy  of  the  slug,  and  from  the  applied  pressure  difference  ( p 2  -  Pj)  and 
the  distance  s  traveled,  the  work  done  (Fig.  17). 

o 

We  use  the  symbol  s^  to  define  the  length  of  the  mixture  column, 
s^  the  length  of  the  final  water  column,  and  s^  the  distance  the  upper 
surface  has  moved.  Then 

s  =  s  -  s,  , 
g  m  f 


s  -as  and  s,  =(1  -  a)  s  /  of. 
g  m  f  ’  g 

The  work  W^  done  in  compression  by  the  driving  air  on  the  water 


column  is 


W1  =  ^2  "  Asg; 


and  the  kinetic  energy  of  the  column  is 

W2  =  (1/2)  Pf  Asf  (Uf)2 


=  A  sf  (p2  -  pj)  a/2,  (1-Q!), 
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where  substitution  was  made  from  Eq.  (11c).  Therefore  E,  the  mechanical 
energy  lost  per  unit  volume  of  vapor  condensed,  is 


E  =  (Wj  -  W2)/  Asg  =  (p2  -  Pj)  £l  -  a  Sf/sg  2  (1  -  or) J  (15) 

=  (P2  -  PL)  /2. 

Thus  one  half  of  the  applied  energy  is  lost  in  collapsing  the  bubbles 
and  heating  the  waiter. 

Of  course,  the  water  is  heated  also  by  the  condensing  steam  giving 
up  its  latent  heat.  (This  heating  accounts  for  the  minimum  pressure  re¬ 
quired  to  collapse  the  bubbles.  )  A  comparison  of  this  heating  effect  with 
the  energy  loss  of  compression  is  instructive. 

The  heat  H  given  up  by  condensing  steam  per  unit  volume  is 
s 

H  =  L  p  , 
s  rg 

where  L  is  the  latent  heat  of  steam  and  p  its  density.  The  heat  H 

rg  e 

produced  by  the  energy  loss  mechanism  is 
He  =  EJ, 

where  J  is  the  mechanical  equivalent  of  heat. 

The  heating  effect  due  to  the  energy  loss  is  thus  a  fraction 

J  (p_  -  p.)/2  p  L  of  the  heat  produced  by  the  condensing  steam.  In  the 

S 

vicinity  of  100° C  it  would  be  necessary  to  have  a  pressure  jump  (p2  -  p^) 
of  about  14  atmospheres  or  200  psi  before  the  two  effects  become  equal. 

The  amplitudes  in  the  investigation  were  confined  to  fractions  of  an 
atmosphere,  (p^  -  p^)  about  10  psi;  therefore  no  attempt  was  made  to  detect 
these  very  small  energy  losses  by  measuring  water  temperature  changes. 
The  heating  due  to  the  energy  loss  mechanism  would  be  the  small  difference 
between  the  total  heating  and  that  of  the  latent  heat  liberation  above.  Initial 
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conditions  were  not  controlled  with  the  precision  that  would  be  required 
for  a  direct  determination  of  the  loss  mechanism  from  the  difference 
between  these  two  relatively  large  quantities. 

III.  EXPERIMENTAL  INVESTIGATION 
A.  Apparatus 

The  propagation  of  ,fshocksM  of  sufficient  amplitude  to  condense 
all  the  vapor  in  a  boiling  water  mixture  was  investigated  with  a  trans¬ 
parent  shock  tube.  The  general  components  of  the  equipment  are  shewn 
diagrammatically  on  Fig.  18  and  pictorially  on  Fig.  19*  The  prominent 
feature  is  a  Pyrex  ’’double  tough”  glass  pipe,  which  is  the  shock  tube 
proper.  It  was  necessary  to  employ  a  glass  pipe  because  visual  obser¬ 
vation  was  needed  to  check  on  uniformity  of  bubble  formation.  The 
pressure  gages  were  mounted  on  annular  brass  inserts  between  glass 
sections*  In  the  illustration,  gage  holders  are  shown  separated  by  6  in. 
of  glass  pipe.  Many  readings  were  taken  with  a  single  insert  having 
a  gage  separation  of  only  2  in. \ Narrow  tubes  are  shown  leading  to  a 
manometer  which, was  used  to  determine  the  bubble  concentration  in 
the  fluid.  The  void  concentration  is  directly  proportional  to  the  height 
difference  in  the  manometer  limbs.  The  difference  in  density  between 
the  manometer  fluid  and  the  water  in  the  pipe  due  to  the  temperature 
difference  was  easily  corrected  by  calibration  with  a  partially  filled 
pipe. 

Bubbles  were  produced  by  boiling  the  water  with  a  wire  heater.  In 
the  reflux  arrangement  shown  a  regular  immersion  heater  boiled  the 
water  near  the  base.  A  separately  adjustable  heater  in  the  region  of 
the  pressure  gages  helped  to  maintain  a  reasonably  uniform  bubble  con¬ 
centration.  The  steam  formed  was  condensed  in  a  water-cooled  condenser 
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and  the  condensate  was  fed  back  to  the  system  near  the  base. 

The  upper  end  of  the  pipe  was  joined  to  the  driver  section.  This  was 
a  steel  pipe  of  the  same  diameter  as  the  main  shock  tube  and  separated 
from  it  by  a  thin  diaphragm.  In  operation,  the  driver  section  was  filled 
with  compressed  air  and  the  diaphragm  was  burst  and  the  shock  wave 
impinged  on  the  surface  of  the  boiling  water. 

The  pressure-time  history  in  the  boiling  water  was  recorded  with 
pressure  gages  on  a  cathode  ray  oscilloscope.  The  traces  of  two  gages 
were  plotted  simultaneously  through  the  use  of  an  electronic  switch 
(Tektronix  53/54C).  The  oscilloscope  time  base  was  started  with  a  third 
transducer  incorporated  in  the  apparatus  very  close  to  the  diaphragm. 

This  assured  the  starting  of  the  trace  well  in  advance  of  the  occurrence 
of  events  in  the  liquid.  Oscilloscope  traces  were  recorded  photographically. 
Details  of  the  pressure  measurement  and  the  diaphragm  are  given  in 
Appendices  F  and  G. 

B.  Photography 

The  events  of  steam  bubble  compression  as  indicated  by  the  oscillo¬ 
scope  traces  occurred  faster  than  the  eye  could  follow.  It  seemed  of 
interest  to  verify  the  sequence  of  events  by  means  of  high-speed  photo¬ 
graphy.  Three  sets  of  pictures  of  these  sequences  are  shewn  in  Fig.  20. 
The  time  scale  shown  below  in  the  first  group  applies  to  all  three.  Differ¬ 
ences  infilm  velocities  have  been  accommodated  by  suitably  mounting  the 
prints.  Water  in  the  column  boiled  at  atmospheric  pressure.  A  shock  wave 
of  12  psi  was  applied  for  the  upper  two  sequences,  and  22  psi  for  the  lower 
one.  We  can  see  very  clearly  on>  these  films  that  as  the  pressure  wave 
travels  down, the  pipe,  steam  bubbles  slowly  collapse.  However,  the  bubble 
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collapse  is  by  no  means  instantaneous,  but  bubbles  are  seen  to  contract 
and  move  down  while  others  further  down  have  already  started  to  collapse. 

In  fact,  a  small  number  of  fairly  large  bubbles  can  be  seen  to  complete 
their  collapse  long  after  those  farther  down  have  been  completely  cleared. 

The  pressure  front  is  difficult  to  identify  on  the  films,  just  as  it 
was  difficult  to  define  on  the  oscilloscope  traces;  definition  is  lost  con¬ 
tinuously  as  the  wave  progresses.  The  velocity  of  the  water  itself  is  quite 
evident,  however,  from  the  bubble  displacement  in  successive  frames. 

These  pictures  confirm  the  analytical  discussions  presented  earlier  in 
the  report, 

C.  Experimental  Results 

On  the  basis  of  the  motion  of  the  liquid- vapor  mixture  under  the 
influence  of  an  incident  shock  wave,  a  sequence  of  pressure  events  has 
been  described  in  Section  II- D.  These  events  were  experimentally  observed, 
including  the  low  propagation  velocity  and  the  pressure  amplification  by 
water  hammer  action. 

The  velocity  of  propagation  of  the  pulses  was  obtained  from  the 
time  difference  of  the  arrival  of  the  pulse  at  the  tw'O  pressure  gages.  It 
has  been  pointed  out  already  that,  in  spite  of  the  rapidly  rising  front  of 
the  shock  wave  striking  the  medium,  the  pressure  rises  slowly  in  the 
mixture.  Moreover,  the  rate  of  rise  decreases  as  the  pulse  progresses 
through  the  fluid.  This  is  brought  out  in  the  upper  oscillogram  of  Fig.  21. 

The  pressure  is  seen  to  rise  at  the  lower  gage  (lower  trace)  about 
0.  8  ms  after  the  start  of  the  pulse  at  the  upper  gage,  whereas  the  upper 
level  or  "knee"  of  the  pulse  is  not  reached  until  2.  2  ms  after  it  was 
reached  at  the  upper  gage.  This  part  of  the  trace  is  seen  to  approach 
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merging  with  the  water-hammer  pulse  at  3.  7  ms.  Complete  merging 

is  sometimes  observed,  as  in  the  lower  pair  of  oscillogram  traces 

O-ik  ^ 

of  Fig.  4,  This  picture  was  taken  with  similar  initial  bubble  concen¬ 
tration,  about  5  per  cent,  but  with  more  amplifier  gain  as  shown  on 
the  ordinate  scale. 

The  "shock"  velocity  is  defined  as  the  distance  between  the  gages 
divided  by  the  time  it  takes  the  front  to  move  from  one  to  the  other. 
Because  different  parts  of  the  front  do  not  have  the  same  travel  time, 
the  concept  of  wave  velocity  is  not  a  precise  one.  However  because  of 
the  very  small  value  of  this  velocity,  even  an  order-of-magnitude  value 
is  of  interest.  Different  velocities  are  of  course  obtained  if  we  measure 
the  velocity  of  the  start  of  the  front  or  the  velocity  of  the  upper  "knee". 
These  different  velocities  have  been  plotted  in  Fig.  22.  The  velocity  of 
the  start  is  shown  as  a  circle,  the  velocity  of  the  "knee"  as  a  triangle 
and  the  midway  points  as  squares.  The  diagram  in  upper  corner  of 
the  graphs  shows  the  pres  sure -time  function  of  the  two  gages;  the  time 
interval  indicated  by  the  appropriate  symbol  divided  into  the  gage 
separation  yields  the  velocity  plotted  on  this  graph  by  the  corresponding 
symbol. 

The  measured  velocity  of  about  50  m/s  (150  ft/s)  at  5  per  cent, 
faling  to  20  m/s  (60  ft/s)  at  50  per  cent  steam  concentration  is  seen  to 
be  somewhat  higher  than  the  computed  velocity  from  Eq.  (11a)  as  shown 
in  the  solid  line. 

Experimental  difficulty  was  predominantly  connected  with  the 
production  of  a  uniform  mixture,  and  measuring  the  bubble  concentration 
at  the*  instant  of  breaking  the  diaphragm.  It  is  seen  that  the  very  low 
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velocity  predicted  is  confirmed,  even  though  it  was  found  difficult  to 
obtain  a  high  degree  of  precision  in  the  measurements. 

The  rate  of  rise  of  pressure  decreases  as  the  pulse  travels  through 
the  medium.  The  rate  of  rise  is,  moreover ,  usually  rather  irregular, 
which  precludes  precise  statement  about  rate  of  the  rise  and  its  change. 

On  the  average  it  was  found  that,  when  the  mixture  was  struck  by  a  true 
shock  on  the  surface,  then' at  a  depth  of  about  1  in.  (2.  5  cm)  the  pressure 
took  about  one  half  millisecond  to  rise  to  its  final  value;  2  in,  (5  cm) 
further  this  rise  time  was  nearly  doubled. 

The  pressure  pulse  traveling  through  boiling  water  finally  produces 
a  large  water  hammer  as  described  in  Section  II- D.  This  is  seen  as  a 
tall  spike  at  3.  7  ms  on  Fig.  21  A.  The  amplitude  of  this  final  pulse  was 
at  first  found  to  be  not  quite  as  large  as  predicted.  A  few  measurements 
were  attempted  subsequently  with  a  longer  column  of  boiling  water. 
Concentration  uniformity  was  poorer.  Fairly  large  concentrations  were 
used  and  averaged  over  the  length  of  the1  column  by  observing  the  increase 
in  height  of  the  column.  Table  3  shows  that  the  long  column  is  needed  to 
observe  the  full  water  hammer  pulse.  The  need  for  the  lcnger  column 
is  readily  understood  when  we  observe  an  expanded  pulse,  Fig.  23.  The 
rise  time  for  the  pulse  is  consistently  about  0.  2  ms,  so  that  when  it  is 
measured  close  to  the  initial  surface,  the  reflection  interferes  destructively 
with  the  initial  pulse.  This  is  very  readily  seen  when  comparing  the 
oscillograms  in  Figs.  23  and  24,  where  a  much  greater  amplifier  gain 
had  to  be  used  for  the  gage  near  the  surface  than  for  the  one  at  the  other 
end  of  the  pipe.  The  correspondence  of  the  upper  trace  of  Fig.  23  with 
Fig.  15A,  and  the* lower  trace  with  Fig.  15C  (inverted  in  Fig.  23  for 
convenience)  confirms  this  behavior. 
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Table  3 


PRESSURE  AMPLIFICATION ,  P/  (p  ,-p  i )  > 

FOR  THE 

WATER  HAMMER 

Concentration.  ^2  ^1 
per  cent  psi 

p/(p2“Pi) 

measured 

P/(P2-Pj) 

computed 

Distance  from 

Surface,  in. 

1 

8 

5 

10 

4 

5 

8 

15 

25 

4 

6 

5 

14 

30 

4 

8 

6 

15 

33 

4 

10 

8 

10 

35 

4 

25 

5 

25 

10 

4 

10 

8 

20 

35 

30 

20 

8 

50 

50 

30 

22 

8 

40 

50 

30 

25 

8 

45 

60 

30 
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IV.  CONSIDERATIONS  OF  FINITE  BUBBLE  COLLAPSE  TIME 


The  experiments  have  led  to  valuable  information  regarding  the 
formation  of  large  amplitude  pressure  pulses  which  can  be  explained 
by  the  formation  of  a  water  hammer.  The  precision  of  experimental 
checks  was  limited  by  the  difficulty  of  obtaining  a  uniform  bubble 
mixture.  The  bubble  concentration  could  be  measured  only  as  an  average 
over  theperiod  of  the  manometer,  whereas  one  really  should  know  the 
instantaneous  value  at  the  time  of  the  shock  wave  passage.  This  instan¬ 
taneous  bubble  concentration  as  plotted  in  Fig.  22  is  not  known  with 
desired  precision.  This  accounts  in  part  for  the  observed  scatter  of 
points  about  the  predicted  curve. 

There  appears  to  be  a  predominant  scattering  above  the  computed 

line.  An  explanation  for  this  phenomenon  was  sought  in  the  finite  bubble 

collapse  times  which  have  been  neglected  in  our  derivation.  The  collapse 

9 

of  a  single  spherical  cavity  was  studied  by  Rayleigh7  using  conservation 

1  2 

of  momentum.  This  was  modified  by  Pies  set  and  Zwick  to  take  into 
account  heat  transfer.  Heat  transfer  was  shown  to  have  only  a  small 
effect  on  the  collapsing  bubble.  This  is  in  contrast  with  the  bubble  growing 
in  a  fluid  superheated  a  few  degrees,  where  the  heat  transfer  slows  up 
bubble  growth  by  an  order  of  magnitude  compared  with  that  predicted  by 
momentum  conservation^’ 

Using  Rayleigh's  approximation  for  the  time  t  required  for  a 
bubble  of  initial  radius  Rq  to  collapse  completely, 

1/2 

to  =  °*  ^15  Ro  C  Pf^p2-pl0 

Substituting  in  this  the  maximum  pressure  in  our  wave  of  10  psi  or  7  x  10 J 


-32- 


Report  No. 
ARF  4132-12 


dyne  cm  ,  a  bubbie  of  0.  5  cm  radius  will  take  about  700  microseconds 
to  collapse*  During  the  collapse  of  the  bubble,  some  pressure  is  trans¬ 
mitted  on  through  the  fluid  to  the  adjacent  region  while  the  bubble  is 
collapsing*  Consequently,  the  next  region  sees  the  pulse  rising  more 
slowly  at  first*  The  next  region  sees  this  initial  rise  further  delayed 
and  so  on  until  the  pulse  takes  a  longer  time  to  rise  to  its  final  value 
than  it  would  take  the  bubble  to  collapse  if  this  pressure  were  suddenly 
applied* 

The  bubbles  in  our  1.  5  in*  pipe  were  less  than  1/ 2  in*  in  diameter, 
usually  much  less,  so  that  a  0,  5  cm  radius  bubble  represents  an  upper 
limit*  The  rise  time  in  our  oscillograms  is  of  the  order  of  1  to  2  milli¬ 
seconds*  The  lack  of  bubble  uniformity  and  consequent  irregular  pressure 
rise  prohibit  a  precise  measurement  of  this  quantity*  Considering  the 
large  order  of  magnitude  by  which  the  velocity  measured  here  differs 
from  typical  velocities  usually  observed  in  single  phase  media,  experiments 
and  theory  are  in  agreement.  The  water  hammer  pulses  had  not  been 
expected,  though  cnce  they  were  observed,  a  very  plausible  explanation 
followed.  The  computation  of  the  amplitude  of  these  pulses  yielded  a 
pressure  amplification  factor  which  was  confirmed  by  measurements  in 
columns  of  adequate  length.  It  must  be  conceded  that  the  precision  of 
these  measurements  was  less  than  that  of  the  earlier  velocity  measure¬ 
ments.  Further  work  on  the  effect  of  bubble  size  would  be  desirable. 

The  generation  of  bubbles  of  uniform  and  controlled  size  is  necessary  for 
the  success  of  such  experiments* 
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V. 


CONCLUSIONS 


Small  amplitude  sound  waves  propagate  through  fluid -vapor  mixtures 
with  very  low  velocities.  The  velocity  depends  on  the  ambient  condition  of 
pressure  (or  temperature)  and  the  relative  quantity  of  the  constituents. 
Further  work  along  these  lines  is  recommended. 

Large  amplitude  waves  in  a  mixture  will  completely  eliminate  one  of 
the  phases  present.  If  the  initial  liquid  content  is  low,  the  liquid  will  be 
evaporated  by  the  heat  of  compression,  with  the  conversion  of  mechanical 
energy  into  heat.  For  larger  fluid  fractions  the  heat  is  absorbed  by  the 
thermal  capacity  of  the  fluid  and  there  is  very  little  temperature  or  pressure 
rise  in  the  mixture;  the  vapor  condenses,  converting  the  potential  energy 
of  the  driving  pressure  partly  into  the  kinetic  energy  of  the  condensed  phase 
and  partly  into  heat,  when  the  inward  velocity  of  the  collapsing  bubble 
surface  vanishes  on  final  collapse.  The  kinetic  energy  of  the  moving  con¬ 
densed  phase  can  give  rise  to  very  high  amplitude  pressure  pulses  upon 
encountering  the  boundary  of  the  low  impedance  mixture. 

On  the  other  hand,  an  explosion  of  finite  energy  content  can  be 
completely  absorbed  in  a  mixture  of  sufficient  volume.  Application  of 
this  phenomenon  is  discussed  in  Section  II-C.  A  system  with  spherical 
geometry  as  distinct  from  the  linear  geometry  considered  in  this  report 
merits  further  investigation. 
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Fig.  4  VELOCITY  OF  SOUND  AS  A  FUNCTION  OF  STEAM  QUALITY 


-39- 


Report  No. 
ARF  4132-12 


VELOCITY  OF  SOUND 


Fig.  5  VELOCITY  OF  SOUND  AS  A  FUNCTION  OF  STEAM  QUALITY 
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)F  SOUND  AS  A  FUNCTION  OF  PRESSURE 


ENTHALPY  1750  8TU/LB  I.I3W.HR/GM 


QNHOS  JO  A1I0013A 


-42- 


Report  No. 
ARF  4132-12 


VELOCITY  OF  SOUNDAT  CONSTANT  ENTHALPY 
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K  Ct=20+3b/2  L 


Fig.  13  WATER  HAMMER  PULSE, 
LONG  MOVING  COLUMN 


Fig.  14  WATER  HAMMER  PULSE, 
SHORT  MOVING  COLUMN 
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Fig.  16  MINIMUM  PRESSURE  OF  PLANE  PRESSURE  FRONT 
REQUIRED  TO  CONDENSE  ALE  THE  STEAM 
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Fig.  17  PRESSURE  PROPAGATION  IN  A  BUBBLE  MIXTURE 


-50- 


Report  No. 
ARF  4132-12 


Fig.  18  DIAGRAM  OF  SHOCK  TUBE 
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Fig.  20  PHOTOGRAPHIC  TIME  SEQUENCE  OF  STEAM  BUBBLE  COLLAPSE 

DURING  A  PRESSURE  TRANSIENT 


Fig.  21  VARIATION  OF  PRESSURE  AT  DISTANCES  1  cm  AND 
6  cm  BELOW  SURFACE;  INITIAL  BUBBLE  CONCEN- 
T RATION,  5  PER  CENT  ' 
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RELATIVE  VOLUME  OF  ST  E  AM,  PE  RCE  NT 

Fig.  22  VELOCITY  OF  PRESSURE  WAVE  IN  BOILING  WATER, 
p  =  10  psi  (0.  7  bars) 
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Fig.  23  ANALYSIS  OF  FINAL  PULSE  NEAR  SURFACE  (a) 
AND  AT  LOWER  END  OF  PIPE~(b) 


Fig.  24  SUCCESSION  OF  PULSES  SHOWING  LONG 
INTERVALS  DUE  TO  FORMATION  BY 
NEGATIVE  PRESSURE 
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Fig.  25  STRETCHING  THE  DIAPHRAGM 
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APPENDIX  A 


DERIVATION  OF  THE  VELOCITY  OF  SOUND  IN  A  WATER- 

STEAM  MIXTURE 


Let 

* 


P  = 

P  = 
V  = 

T  = 
S  = 

q  = 

Subscripts:  f 

g 

e 


velocity  of  sound, 

pressure, 

density, 

specific  volume  =  1/p, 

temperature, 

entropy, 

steam  quality, 

mass  of  steam/total  mass, 

refers  to  the  fluid  (water) 

refers  to  the  gas  (steam) 

refers  to  the  change  of  parameter  on  evaporation. 


Thus 


The  basic  equations  are 

2  ,dP\  v2//aV\ 

c  =  (*£>  =  ~  V  /<^>  > 


dS  =  ^  dp  +  ^  dq 


dV  =  (S  dp  +  &  dq 

p  q  4  p 


(A-l) 

(A- 2) 

(A-  3) 


S  =  q  Sg  +  (1-q)  Sf  =  q  Se  +  Sf 


(A-4) 


.A-l 
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V  =  q  v  +  (1-q)  vf  =  ,  Ve  +  V£  . 


(A- 5) 


Since  in  (A-l)  S  =  constant,  dS  =  0.  Substituting  this  in  (A-2)  and 


rearranging. 


d9  -  "  (■grj) 

(dS=0)  P  q 


(A- 6) 


Now  differentiate  (A-4)  and  (A-5),  bearing  in  mind  that  V  ,  V ,,  s. . 

6  I 

Se>  and  Sf  are  not  functions  of  q  but  are  uniquely  related  to  pressure: 


_0  ,  ,9S,  _  dSe  ^  dSf 

Wp  "  Se’  and  ^  ~  q  Sp'  +  3p~  ; 


(A-7) 


3  y  3  y  dVe  dVf 

(^q)p  ■  ve’  and  (^p}q  ■  q  dF^  +  dF 


(A-8) 


Substitute  (A-4)  and  (A-7)  in  (A-6);  then 


dq  =  -  (1/S  )  q  dS  /dp  +  dS,/dp  dp. 
(dS=0)  e  e  1 

Substitute  (A-8)  and  (A-9)  in  (A- 3);  then 


(A-9) 


dV  dV,  V  dS  dS,  ] 

i  *  *  e  .  1  6/  C  ,  X\i  j 

dV  =  9  -^r-  +  -3TT  -  c—  (q  -x=-  +  -T3— ) !  dp. 

(dS=0)  dp  dp  Se  dp  dp  i 


Collect  the  terms,  so  that 


3V  dV  V  dS  dV  V  dS 

Vs  =  q  (-dF'  "  s;  -7IF  }  +  “FT  ■  ST  "dF  * 


(A  -  10) 


( A  -  1 1 ) 


then 


,2  _  ,9p 


Vs 


(A- 12) 
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that  is 


Table  A  -  1 


IMPEDANCE  OF  BOILING  WATER 


p 

psi 

Impedance 

,,  ,*2  -1  -2  -1 
lb  ft  sec  gm  cm  sec 

P 

bars 

1 

20.  5 

10.0 

0.  068 

10 

158 

77.  8 

.  68 

20 

294 

142 

1.  36 

50 

625 

308 

3.4 

100 

1,050 

517 

6.  8 

200 

1,960 

965 

13.  6 

500 

4,070 

2,000 

34 

1000 

6,800 

3,  340 

68 

2000 

11,400 

5,610 

136 

3000 

13, 300 

6,  540 

340 

3206 

—  20,  000 

-10,000 

364 

A-4 
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We  may  simplify  furthermore  by  considering  that 


then 


=  (Ve/Se)(dSf/dT)(dT/dp) 
=  <Ve/Se)(Cf/T)(Ve/Se), 


(A- 15) 


(A- 16) 


where  we  have  defined  the  specific  heat  of  the  fluid  =  dQ/Dt  .  = 

TdS^/dT,  and  have  made  use  of  the  Clausius -Clapeyr on  relation  dT/dp  = 

V  /S  .  If  we  further  define  a  latent  heat  L  =  TS  ,  then 
e  e  e 


pc 


=  Se/(Ve  V  Cf/T) 

=  L/(VeVTCf)  • 


(A- 17) 


The  flow  rate  through  a  pipe  is  limited  to  the  sonic  velocity.  The 
rate  of  mass  flow  of  fluid  under  these  so-called  "choked"  conditions  is 
Ape  =  AL/(Ve*V  TCf),  whew-  A  is  the  area  of  the  pipe. 


A-5 
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APPENDIX  B 


USE  OF  STEAM  TABLES  TO  COMPUTE  THE  VELOCITY  OF  SOUND 


The  quantities  V  ,  Vf,  Se  and  Sf  required  by  the  formula  derived 

3 

in  Appendix  A  are  found  tabulated  in  Keenan  and  Keyes'  steam  tables  . 
The  derivative  of  the  quantities  is  obtained  by  a  linear  approximation 


from  adjacent  lines  in  the  table,  thus; 


d.V  (V  )  . 

!  e\  ~  e'n  +  a. 

'dp  ;  “  p  . 

r  n  n  +  a 


(V  )  . 

'  e  n-b 


(B-l) 


where  p  ,  etc.  ,  refer  to  the  values  in  the  nth  lines  of  the  table.  In 
rn 

general,  a  and  b  were  taken  as  unity;  that  is,  values  were  computed  from 
adjacent  lines  in  the  table  except  where  the  intervals  were  suddenly 
changed.  In  that  case,  a  or  b  was  chosen  so  that 

Pn  +  a  -  pn  =  pn  -  Pn-b  •  (B'2) 


The  error  introduced  by  this  technique  is  partly  due  to  the  fact 

tha,t  the  slope  of  the  chord  differs  slightly  from  the  slope  of  the  tangent 

at  the  mid-point  of  the  arc,  and  partly  due’ to  the  small  differences  between 

adjacent  variables*  Curve  fitting  with  approximate  analytical  expressions 

and  differentiating  eliminates  the  first  error,  but  not  the  second.  That 

is,  we  cannot  fit  any  curve  any  more  accurately  than  the  given  points. 

The  errors  due  to  the  difference  between  chord  and  tangent  can  be  shown 

2 

to  be  of  the  order  of  (Ap/p)  for  the  functions  of  interest  (see  Appendix  C). 

2 

Table  B-l  gives  a  few  values  of  (Ap/p)  and,  the  percentage  differencVof 
the  last  digits  of  A  V  ,  AS^,  anc*  A  *  The  precision  of  these  last 
variables  is  the  rounding- off  in  the  last  digit.  In  retrospect  it  appears 


B-l 
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that  the  precision  could  have  been  increased  somewhat  by  lengthening  the 
chords,  thus  increasing  the  values  of  AV  ,  AS  ,  and  AS,.  and  thereby 
decreasing  the  uncertainty  in  the  last  digit.  This  is  especially  true  in 
the  region  of  100  psi.  The  optimum  length  of  chord  would  be  that  for 
which  the  two  errors  here  discussed  are  approximately  equal.  It  was 
not  considered  worthwhile  to  repeat  the  calculations  for  this  minor 
improvement  of  precision. 

Table  B-l 

ERROR  DUE  TO  CHORD  APPROXIMATION 


P  -2 
lb  in. 

Ap 

( Ap)2 
p 

per  cent 

AV 

e 

per  cent 

AS 

e 

per  cent 

ASf 

per  cent 

1 

0.  2 

4.  0 

0.  07 

0,  2 

0.  5 

2 

0.  2 

1.  0 

0.  03 

0.  3 

0.  8 

5 

0.  5 

1.  0 

0.  10 

0.  3 

0.  7 

10 

1.  0 

1.  0 

0.  15 

0.  3 

0.  6 

20 

1.  0 

0.  25 

0.  01 

0.  7 

1.  4 

50 

1.  0 

0.  04 

0.  3 

1.  7 

3.  0 

100 

1.  0 

0.  01 

1.  2 

3.  0 

5.  0 

200 

10.  0 

0.  25 

0.  9 

1.  0 

2.  0 

500 

20.  0 

0.  16 

0.  1 

0.  6 

1.  3 

1000 

100.  0 

1.  00 

0.  1 

0.  2 

0.  3 

2000 

100.  0 

0.  25 

0.  5 

0.  2 

0.5 
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APPENDIX  C 


ERRORS  INTRODUCED  BY  USING  THE  CHORD  APPROXI- 
MATION  FROM  STEAM  TABLES 


First  consider  the  difference  of  slope  between  the  chord  and  the 
tangent.  The  chord  is  chosen  to  be  symmetrical  about  the  point  of  interest 
as  defined  in  Eq.  (B-2).  Over  a  limited  region  the  value  of  Ve  varies 
approximately  as  1/  p  and  Sg  varies  approximately  as  log  p.  For  both 

2 

these  types  of  function  we  shall  now  show  the  error  to  be  equal  to  (Ap/  p)  , 


*  /  a  n  f  t  nAx 

y- Ay  =  (x-Ax)  =  x  1  -  — —  +  -  - 

n  F  0  nAx  0  n(n-l)(n-2)  ,Axx3 
2Ay  =  x  ^2  +  2  2  .  g  (— )  + 

=  [1+  (^)2  +  -- 


Now  the  definition  of  the  tangent  is  the  slope  of  the  chord  as  Ax  — *-  0; 

2 

thus  (1  +  ^n"  (— )  +  higher  powers  of  ^ )  is  the  error  in  using 

a  chord  symmetrically  disposed  about  the  point  of  interest  Now  V  varies 
approximately  as  p Hence,  the  chord  symmetrically  disposed  about  the 
point  of  contact  differs  from  the  tangent  at  that  point  by  factor  differing 
from  unity  by  • 

Case  2  , 

y  =  log  x, 

Ay  =  log  (x  +  Ax)  -  log  x, 

Ay  =  log  (1  +  A x/x), 
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-Ay  =  log  (1  -  Ax/  x); 
therefore 

2Ay  =  log  (1  +  A x/x)  -  log  (1  -  Ax/ x) 

=  Ax/x  +  (Ax/x)2  +  (Ax/x)3  +  —  ~  +  Ax/  x  -  (Ax/x)2  +  (Ax/x)3  + 

=  2  Ax/  x  +  2  (Ax/  x)3  +  ---; 

hence 


Now  Sg  and  approximate  a  log  p  function  over  a  limited  regions,  so  here 

2 

again  chord  errors  are  of  the  form  (Ap/  p)  . 

It  should  be  noted  that  if  a  non- symmetrical  chord  is  chosen  to 
approximate  the  tangent,  so  that  the  one  end  of  the  chord  is  at  a  point 
of  interest,  then  the  error  would  have  been  Ap/p  instead  of  (Ap/  p)  , 

This  is  seen  immediately  in  the  above  derivation,  since  the  first  terms 
in  the  binomial  expansion  cancel  owing  to  the  symmetrical  disposition  of 
the  chord. 


C-2 
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APPENDIX  D 


THE  RELATIONSHIP  BETWEEN  VOLUME  AND  PRESSURE 
ALONG  A  LINE  OF  CONSTANT  ENTROPY 


The  expression  for  volume  in  terms  of  pressure  (Section  IV)  may¬ 
be  obtained  as  follows: 

As  in  Appendix  A,  we  write 


v  =  q  Ve  +  Vj  , 

q  =  (V-  Vf)/Ve, 

(D-  1) 

S  =  q  S,  +  S  . 

q  =  (S  -  Sf)/Se  ; 

(D-2) 

thus 

V  =  (S  -  sf)  Ve/se  +  vf , 

=  S  V  /S  -  (Sr  vjsa  - 
e  e  i  e  e 

Vf), 

=  SX^Xr 

(D-  3) 

The  parameters  Vq/S^  =  X  ,  (p) 

(D-4) 

and  Sf  Ve/Se  -  Vf  =  X£  (p) 

(D-5) 

were  plotted  as  a  function  of  pressure. 

Analytical  expressions  for  X2  are  quite  awkward,  but  by  changing 
the  arbitrary  zero-reference  for  the  specific  entropy  and  calling  this 
modified  function  X^,  simple  power  laws 


X  L  (p)  =  Ap "  a  +  Bp”b, 

(D-6) 

X3  (p)  =  Fp~f  , 

(D-7) 

can  be  fitted  to  the  functions 

X1  <P>  *  Ve/Se  > 

and 

x3  (P)  =  (Ve/Se)  (Sf  +  SQ)  -  Vf  . 

(D-8) 

D-  1  Report  No, 

ARF  4132-12 


Substituting  Eq.  (D-8)  in  (D-3),  we  obtain 


v  =  (s  +  so)  ve/se  -  (Sf  +  so)  ve  se  -  Vf 


=  (Ap"a  +  Bp"*5)  (S  +  S  )  -  Fp“f. 


For  the  units  in  available  tables,  pressure  in  lb  in.  ~2  and  specific  volume 


3  - 1 

in  ft  lb  ,  the  constants  best  fitting  the  curves  were 


A  =  113, 


a  =  1. 


=  0.  24. 
o 


B  =  67.  5, 
b  =  0.  7, 


F  =  67.  5  , 
f  =  0.  7, 


The  equality  of  B  and  F  and  of  b  and  f  then  permitted  further  simplifi- 
cation  to  the  final  form, 


V  =  Aq  (S  +  So)  p-1  +Bo(S  +  So-  1)  p- 


(D-9) 


Where  A  -  Aq  and  B  -  Bq.  Numerical  values  of  the  constants  in  Eq.  (D-9) 
are  given  in  Table  2  for  various  systems  of  units. 
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APPENDIX  E 


MINIMUM  PRESSURE  NEEDED  TO  CONDENSE  ALL  THE  STEAM 


To  obtain  the  minimum  pressure  needed  to  condense  all  the  steam 
in  a  boiling  mixture,  we  approach  the  problem  from  the  other  end  by 
considering  the  small  amplitude  region  and  letting  the  pressure  amplitude 
increase  until  we  reach  the  point  of  complete  condensation.  In  addition, 
we  make  use  of  the  observation  in  Section  II-C  that  the  Rankine-Hugoniot 
in  the  boiling  mixture  is  very  close  to  the  adiabat. 

Next  define  the  volume  concentration  of  steam,  o?  ,  and  express  it 
in  terms  of  the  mass  concentration,  q: 


(1/0!  -  1)  Vg  =  (l/q-  l)Vf. 

(E-l) 

For  boiling  water  a  <  0,  9;  mass  concentration  of  steam,  q, 

so  that  we  may  write  approximately 

is  very  small 

q  S  o;Vf/  [(1  -  0!  )  Vg], 

(E-2) 

and 

d  q/dffi  =  Vf/Jjl-0!)2  Vg]  . 

(E-3) 

Now  for  q  «  1,  Eq.  (A-9)  becomes 

Op/3q)s  S  -  Se/(dSf/dp). 

(E-4) 

Combining  (E-3)  and  (E-4), 

(8p/  9  0!)s  =  -  Vf  Se/  [vg  (1  -O')2  (dSf/dp)]  . 

(E-5). 

E-l 
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To  obtain  the  maximum  pressure  which  just  barely  does  not  condense  all 
the  vapor,  we  integrate  from  the  existing  concentration  &  to  zero  and 

obtain:  o 

/> 


Pmax  *  -  (  ( 1  -  a)' ' 2  Vf  Se  da  /  [Vg  <dS£/  dp)] 

l 

(E-6) 

=  Of  (l  -  a)" 1  vf  Se/  [ vg  (dsf  /dp)  ]  . 

(E-7) 

This  is  now  of  course  also  equal  to  the  minimum  "shock"  pressure 

necessary  to  condense  all  vapor.  Its  value  is  plotted  in  Fig.  E-l.  Now 

let  us  substitute  this  minimum  pressure  into  the  Eq.  (11a)  to  obtain  the 

velocity  of  a  pulse  barely  strong  enough  to  condense  all  the  vapor. 

<Ul2)min  =  <?2  "  Pl>min  V  O  (1  “ 

(E-8) 

=  (l  -  a)"2  vf2  se  /  [vg  (dSf/  dp)] 

(E-9) 

=  V2  Se/[Vg  (dSf/dp)J. 

(E-  10) 

This  is  very  similar  in  magnitude  to  the  velocity  of  sound  in  the  mixture, 
Eq.  (A-13),  using  the  following  relations : 

V  =  q  Ve  +  Vf  =Vf/(l  -  oi  ), 

and 

q  d  Ve/dp  -  q(Ve/Se)(dSe/dp)  +  d  Vf/dp«  ( Ve/Se)(dSf/dp)  . 

Thus  the  velocity  of  propagation  of  pulses  of  amplitude  just  sufficient  to 
condense  al'l.the  steam  is  approximately  the  sonic  velocity. 


E-2 
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APPENDIX  F 


PRESSURE  MEASUREMENTS 

The  pressure  was  measured  with  quartz  transducers  (Kistler  Instru¬ 
ment  Company,  Tonawanda,  New  York,  Type  601).  These  gages  have  an 
overall  diameter  of  0.  25  in.  and  an  active  diaphragm  of  about  0.  12  in. 

The  self-resonance  frequency  is  250,000  cps.  A  check  with  a  shock 
wave  in  air  traveling  across  the  diaphragm  showed  that  the  time  for  the 
output  to  reach  90  per  cent  of  its  final  value  was  less  than  5  microseconds. 

The  gages  were  calibrated  by  two  independent  techniques.  An  electro 
meter  calibrator,  Kistler  PZ-6,  was  used  to  measure  the  charge  produced 
on  the  crystal  for  a  static  pressure  produced  by  a  dead  weight  tester. 

An  acoustic  technique  was  also  used.  In  this,  the  gage  was  subjected 
to  a  1  kc  alternating  pressure  produced  by  a  loudspeaker  in  an  air-filled 
cavity.  The  pressure  was  measured  with  a  condenser  microphone,  for 
which  a  reciprocity  calibration  was  available. 

The  pressure  used  in  the  first  calibration  was  500  psi;  that  used 
in  the  second  was  0.  15  psi.  The  sensitivity  achieved  with  the  former 
technique  was  about  10  per  cent  greater  than  that  obtained  with  the  latter. 
The  precision  of  the  latter  technique  is  not  expected  to  be  better  than  that 
achieved;  however,  it  provides  a  good  check  for  the  operation  of  the  gage 
at  audio  frequency. 

It  was  found  convenient  to  connect  the  gages  via  cathode  followers 
to  the  separate  inputs  of  a  switching  amplifier  (Tektronix  Type  53/54C) 
to  view  the  signals  from  the  two  stations  separately.  Special  short 
cables  (about  4  in.  )  were  used,  so  that  the  load  on  the  crystals  was 
minimized  (20  pp  f ,  700  megohms).  In  this  way  the  voltage  sensitivity 
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(0.  1  v/cm)  of  the  oscilloscope  was  adequate  without  further  amplification. 
An  earlier  scheme  using  the  output  of  both  gages  on  the  same  channel 
worked  well  in  a  single  phase  medium  (air),  but  led  to  unsatisfactory- 
pictures  for  pulses  traveling  through  the  mixture. 
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APPENDIX  G 


THE  DIAPHRAGM 


Some  time  had  to  be  devoted  to  a  search  for  a  suitable  diaphragm 
material,  since  standard  plastic  materials  do  not  perform  well  at  the 
temperature  of  boiling  water.  On  the  other  hand  operation  at  lower 
temperature  and  greatly  reduced  pressure  proved  to  increase  the 
difficulties  in  producing  uniform  mixtures.  Tests  with  neoprene  0.  015 
in.  thick  and  0.  031  in.  thick  were  quite  satisfactory.  The  material  was 
stretched  as  far  and  as  uniformly  as  possible,  as  shown  in  Fig.  25,  on 
a  ring  slipped  over  a  jig.  A  standard  hose  clamp  held  the  material  in 
place  (Fig.  26).  By  the  pre stretching  technique  used,  rapid  opening  of 
the  diaphragm  to  well  over  three  fourths  of  the  ring  area  was  assured. 

By  making  the  ring  2  in.  in  diameter  for  use  in  conjunction  with  the 
1-1/2  in.  I.  D.  pipe,  the  remaining  parts  of  the  diaphragm  presented 
no  severe  obstruction.  A  slight  area  change  near  the  diaphragm  was 
also  required  by  the  steam  outlet.  This,  it  was  later  shown,  could  have 
been  greatly  reduced.  Small  oscillations  due  to  turbulence  in  this  region 
were  not  very  important  because  the  time  scale  was  so  long  that  these 
turbulence -induced  fluctuations  were  virtually  completely  damped  out. 

A  few  observations  of  shock  waves  in  air  at  a  point  14  in.  below  the 
diaphragm  indicated  a  well-formed  shock  wave  at  that  point  (rise  time, 
as  well  as  could  be  observed  with  available  instruments,  appeared  to 
be  less  than  5  microseconds).  The  strength  of  the  neoprene  at  elevated 
temperatures  decreased  with  time,  and  if  heated  for  a  sufficiently  long 
time,  would  sometimes  burst  without  any  applied  air  pressure.  However, 
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for  a  few  minutes  the  diaphragm  would  withstand  over  75  per  cent  of 
its  room-temperature  bursting  pressure  (12  to  16  psi  for  the  thin 
material,  and  2  2to  3  0  psi  for  the  thicker  one). 

The  use  of  neoprene  rather  than  cellulose  acetate  or  mylar  as 
a  diaphragm  material  has  the  advantage  that  the  diaphragm  ruptures 
cleanly  and  does  not  shatter,  thereby  maintaining  the  shock  tube  and 
valve  system  free  of  fragments  and  reducing  the  cleaning  problem. 
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APPENDIX  H 


THE  EQUILIBRIUM  ASSUMPTION 

The  assumption  of  equilibrium  conditions  made  at  the  outset  of  the 
analysis  (Section  II-A)  is  valid  for  sufficiently  low  frequencies  and  for 
particles  of  the  discontinuous  phase  (that  is,  water  droplets  in  steam, 
or  steam  bubbles  in  water)  of  sufficiently  small  size.  From  thermal 
diffusion  we  may  estimate  the  frequency  limit  for  particles  of  a  given 
size,  or  the  particle  size  limit  for  a  given  frequency,  so  that  equilibrium 
conditions  still  hold. 

15 

Jakob  and  Hawkins:  Elements  of  Heat  Transfer  gives  in  Fig.  IV- 7 
the  temperature  at  the  center: of;  a  sphere  whose  surface  temperature  is 
suddenly  raised.  This  curve  plots  (T  -  T  )/(T  -  T  ),  the  ratio  of 

CL  S  CO  s 

the  difference  between  the  center  temperature  T^  at  time  t  and  the 

surface  temperature  Tg,  to  the  initial  difference  in  temperature  (Tcq  -  Tg 

at  time  t  =  0,  when  the  surface  temperature  was  raised  from  Tco  to  Tg  . 

A  few  points  taken  from  this  graph  have  been  tabulated  in  Table  H-l  in 

2 

terms  of  the  dimensionless  quantity  4ctt/D  where  a  is  the  thermal 
diffusivity  of  the  sphere  and  D  is  its  diameter. 


Table  H-l 

CENTER  TEMPERATURE  OF  A  SPHERE 


4  a.  t/ 

0.  1 

0.  15 

o 

o 

. 

0.  3 

0.  4 

0.  5 

T  .  -  T 
ct  s 

T  -  T 

CO  s 

0.  7 

0.4 

0.  25 

0.  1 

0.  04 

0.  01 

H-l 


Report  No. 
ARF  4132-12 


From  this  table  we  see  that  most  of  the  heat  transfer  has  taken  place 
in  a  time  such  that  4  crt/D2  has  values  between  0.  1  and  0.  5.  Without 
going  into  the  details  of  the  temperature  distribution  throughout  the  sphere 
at  each  instant  and  calculating  the  heat  fluxes,  we  see  that  a  very  large 
fraction  of  the  total  heat  transfer  has  taken  place  when  4  art/D  has  the 
value  0,  2.  From  this  we  obtain  a  time  constant  t  =  .  05  D  /a  in  which 
we  may  say  all  heat  transfer  has  taken  place.  The  reciprocal  of  this 
time  constant  will  give  us  the  value  of  the  frequency  below  which  the 
equilibrium  assumption  is  satisfactory.  This  frequency  has  been  tabulated 
in  Table  H-2  for  liquid  drops  in  vapor  and  vapor  bubbles  in  liquid. 

Table  H-2 

FREQUENCY  LIMIT  FOR  EQUILIBRIUM  ASSUMPTION 


Diameter, 

Fre  quency  Limit 

For 

mm:  1*  0 

0.  1 

Frequency, 

0.01 

cps 

water  at  32° 

* 

00 

280 

28,000 

water  at  400° 

3.4 

340 

34,000 

saturated  vapor  at  32° 

31,000 

3.  1  x  106 

3.  1  x  108 

/ 

saturated  vapor  at  200* 

340 

3.  4  x  104 

3.  4  x  10° 

saturated  vapor  at  400° 

25 

2.  5  x  103 

2.  5  x  105 

A 

saturated  vapor  at  600° 

5.  5 

5.  5  x  102 

5.  5  x  10 
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